NPS-EC-93-022 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


94-09402 

iiiiinii 


THESIS 


DTI. 

1  elects 


Low  Altitude  Near-the-Horizon  Propagation: 

A  Comparison  Between  RPO  and  M-Layer 

by 

Chi-Wei  Wu 
December  1993 

Approved  for  public  release;  distribution  is  unlimited. 

Prepared  for:  Program  Executive  Office,  Theater  Air  Defense  (D2) 
2531  Jefferson  Davis  Highway 
Arlington,  VA  22242-5170 

94  3  25  084 


f 


one  QUAUr^  il^&PEClEO  1 


NAVAL  POSTGRADUATE  SCHOOL 

Monterry,  California,  93943 


Rear  Admiral  T.A.  Mercer  H.  Shull 

Superintendent  Provost 


This  thesis  was  prepared  in  conjunction  with  research  sponsored  in  part 
by  Program  Executive  Office,  Theater  Air  Defense  (D2)  under  NPS-EC-93>022. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 


Released  by: 


P.J.  Mi 

''1 

Dean  of  1 

Research 

REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704^188 

AiMe  (•poAngburMn  tor  MM  oeNwttcn  of  nl«mMiDni«MlmaM  1  hour  par  iMponM.  ndudngtw  ton*  tor  (•VMonamMrMtarto.bMMhnaaKCnadatoMuraM 

OOtwring  «nd  ntontoomo  tto  Ohw  noopop.  ond  ccmptoMno  wd  towuuwng  tm  cottaon  N  oitomtoian.  Sand  oommoms  ipgaiding  Ihiu  burdun  oMimto  or^  oifw  awM  at  tip 

*y*r5-^  .  Owogitoto  lor  I^matotn  Oparbiionb  and  RaparwjgtSjStoiaon 

Obwa  HigtMay.  Suaa  1204.  ArlngAn,  VA  22202-4102.  and  to  (to  Ottioa  o(  Mattogarnaru  and  Budgol.  ^pantorb  Baductton  Pto|act(0704%TM),  Wtotoupp^  Ob  20(01 

1.  AQENCY  USE  ONLY  (Laava  ManA)  I  2.REPOnT  DATE  j.  REPORT  TYPE  AND  OATES  COVERED 

1  December  1993  Technical  Report/Master’s  Thesis 

4.  TITLE  AND  SUBTTFLE 

Low  Altitude  Near  the  Horizon  Propagation: 

A  Comparison  Between  RPO  and  M-uyer 

S.  FUNDING  NUMBERS 

6.  AUTHOn(S) 

Chi-Wei  Wu 

7.  PERFONMINO  ORGANIZATION  NAME(S)  ANC;  A00RESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

9.  PERFORMING  0RGAMZAT10N 
REPORT  NUMBER 

NPS-EC-93-022 

9.  SPONSORINGMOMTORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Program  Executive  Office,  Theater  Air  Defense  (D-21) 

Attn:  Mr.  George  Hamilton 

2531  Jefferson  Davis  Highway 

Arlington.  VA  22242-5170 

10.  SPONSORING/MONITORING 

M&tCi  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTE 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or 
position  of  the  Department  of  Defense  or  the  United  States  Govemment 

12a.  DISTRIBUnONIAVAILABIUTY  STATEMENT 

Approved  for  public  release; 

Distribution  is  unlimited. 

129.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Predictions  of  propagation  loss  made  by  the  computer  programs  Radio  Physical  Optics  (RPO) 
Computer  Software  Configuration  Item  (CSCI)  and  M-Layer  are  compared.  The  results  of  the 
high  frequency  parabolic  equation  approximatioivas  formulated  in  RK^  agree  almost  always 
with  those  derived  from  the  low  frequency  modal  computation  as  formulated  in  M-Layer.  But 
at  low  altitudes  in  the  neighborhood  of  the  radar  horizon,  deviations  between  RPO  and  M-Lay¬ 
er  become  significant  for  some  cases.  RPO  appears  not  to  be  able  to  prop>erly  account  for  the 
fects  of  a  high  altitude  surfiice-based  duct  at  a  short  range.  Since  the  discrepancies  fall  in  regions 
of  importance  to  naval  operations,  a  definitive  resolution  is  an  uig^t  task  to  be  undertaken  in 
the  immediate  future. 

14.  SUBJECT  TERMS 

Propagetion  loss,  Radi( 
Waveguide  Mode,  Eva] 

)  Physical  (Dptics  CSCI,  M-Layer,  Parabolic  Equation, 
loration  Duct,  Surface-Based  Duct 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 

19.  SECURTTY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 
OFTMSPAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED 

2D.  LMrrATION  OF  ABSTRACT 

UL 

NSN  754001-280-Sa00 


1 


Standard  Fbim  298  <R*v.  2-89) 
PTMOttad  br  ANSI  SU  Z39-1B 


Approved  for  public  release;  distribution  is  unlimited 


Author 


Ai^Hoved  by: 


LOW  ALTITUDE  NEAR  THE  HORIZON  PROPAGATION: 
A  COMPARISON  BETWEEN  RPO  AND  M-LAYER 

by 

Chi-Wd  Wu 

Lt,  Republic  of  China  Navy 
Chung  Cheng  Institute  of  Technology,  1987 

Submitted  in  partial  fulfillment  of  the 
lequiiemoit  for  foe  degree  of 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  foe 

NAVAL  POSTGRADUATE  SCHOOL 
December  1993 

Chi'Wd  Wu 


Hung'Mou  Lee,  Thesis  Advisor 


Michad  A.  Morgan,  Chairman 
Department  of  Electrical  and  Computer  Engineering 


ii 


ABSTRACT 


Predictions  of  propagation  loss  made  by  the  romputer  programs  Radio  Physical  Optics  (RPO) 
Computer  Software  Configuration  Item  (CSCI)  and  M-Layer  are  compared.  The  results  of  the  high 
frequency  parabolic  equation  approximation,  as  formulated  in  RPO.  agree  almost  always  with 
those  derived  from  the  low  frequency  modal  computation  as  formulated  in  M-Layer.  But  at  low 
altitudes  in  the  nei^bothood  of  the  radar  horizon,  deviations  between  RPO  and  M-Layer  become 
significant  for  some  cases.  RPO  appears  not  to  be  able  to  properly  account  for  the  effects  of  a 
high  altitude  surface-based  duct  at  a  short  range.  Since  the  discrepancies  fall  in  regions  of 
importarKe  to  naval  operations,  a  definitive  resolutitm  is  an  urgent  ta^  to  be  undertaken  in  the 
immediate  future. 
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I.  INTRODUCTION 


A.  THEORETICAL  BACKGROUND 

The  importance  of  environmental  effects  on  communication  links  and  radar  systems 
has  been  recognized  for  well  over  two  decades.  The  proliferation  of  computer  programs 
such  as  the  Integrated  Refraction  Effects  Prediction  System  (IREPS)  [Ref.  1]  and  the 
Engineer’s  Refractive  Effects  Prediction  System  (EREPS)  [Ref.  2]  within  the  Navy  for 
the  prediction  of  propagation  loss  of  radio  waves,  together  witii  the  general  availability 
of  computing  power  brought  along  by  the  personal  computer  (PC)  revolution,  has  greatly 
increased  the  Navy’s  awareness  of  such  effects. 

EREPS  essentially  is  a  PC  version  of  IREPS,  with  provisions  for  greater  flexibility 
in  setting  input  parameters.  It  employs  a  combination  of  different  methodologies  to 
deduce  the  propagation  factor,  ray-optics  within  line-of-sight;  curve-fitting,  together  with 
frequency  scaling  based  on  the  assumption  that  a  single  mode  contributes  to  the  complete 
field  strength  [Ref.  3],  output  from  M-Layer  [Ref.  4,  5]  in  the  over-the-horizon  region, 
and  linear  interpolation  of  the  results  from  ray-optics  and  M-Layer  in  the  region  in 
between.  Through  such  curve-Etting  techniques,  EREPS  (and  IREPS)  gains  speed  at  the 
expense  of  accuracy. 

The  M-Layer  program,  though  time-consuming  to  run,  provides  results  beyond  the 
crude  approximation  of  EREPS,  especially  in  the  penumbra  region  where  hundreds  of 
modes  are  required  to  give  an  accurate  reading.  As  can  be  seen  in  Chapter  m,  M-Layer 
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works  well  at  low  altitudes,  even  in  the  illuminated  region.  The  theory  of  M-Layer  is 
formulated  in  terms  of  the  excitation  and  propagation  of  electromagnetic  waves  in  the 
earth-atmospheric  waveguide.  The  earth,  in  fact,  is  treated  as  a  flat,  lossy,  homogeneous 
half-space.  Its  curvature  is  compensated  for  by  postulating  the  existence  of  a  linearly 
increasing  component  of  the  index  of  refraction  in  the  atmosphere,  in  addition  to  the 
slight  natural  variation  of  the  index  of  refraction  of  real  air.  The  rate  of  increase  of  the 
linear  component  is  set  to  the  inverse  of  the  effective  earth  radius.  This  new  index  of 
refraction  is  called  the  modified  index  of  refi:action,  and  this  method  of  simplifying  the 
problem  constitutes  the  earth-flattening  approximation  [Ref.  6].  Combined  with  the 
piecewise-linear  approximation  to  the  modified  index  of  refraction,  the  eigenfunctions  of 
the  flattened  earth-atmospheric  waveguide  are  Airy  functions  of  different  variables  over 
different  regions  of  the  atmosphere  of  linearly  varying  modified  refractive  indices.  The 
waveguide  modes  are  found  numerically,  which  is  the  process  taking  up  most  of  the 
computation  time. 

Recently,  the  RDT&E  Division  of  the  Naval  Command,  Control  and  Ocean 
Surveillance  Center  (NRaD),  which  previously  produced  IREPS,  EREPS  and  M-Layer 
when  it  was  known  as  Naval  Ocean  Systems  Center  (NOSQ.  came  out  with  a  Radio 
Physical  Optic  (RPO)  program  [Ref.  7].  It  is  structurul  in  the  same  spirit  as  EREPS: 
different  approximations  are  used  over  difforent  regions  of  space.  The  pan  in  EREPS 
which  relies  on  curve-fitting  of  M-Layer  predictions  is  replaced  with  direct  computation 
based  on  the  parabolic  equation  approximation  (PE)  (Ref.  8,  9].  This  makes  RPO  far 
superior  to  EREPS.  Computational  expediency  precludes  extensive  application  of  PE. 
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Instead,  PE  computations  are  launched  at  a  range  of  2.5  km.  unless  the  incident  ray  to  the 
earth  gets  near  the  diffraction  limit  first.  A  split-step,  outward  advancing  scheme  for 
solving  the  parabolic  equation  is  utilized.  At  each  step,  a  fast  Fourier  transform  (FFT)  is 
performed.  The  FFT  size  is  set  between  7  to  10  powers  of  2.  while  the  separation  between 
FFT  data  points  is  limited  to  within  2  to  49  wavelengths.  Hence,  PE  computations  will 
never  be  done  above  a  height  of  50,176  wavelengths.  At  10  GHz,  this  represents  a  height 
of  no  more  than  1.50528  km.  To  avoid  aliasing,  the  FFT  coefficients  computed  from  3/4 
of  the  height  to  the  top  of  the  region  where  the  FFT  is  carried  out  are  filtered  with  a 
factor  which  decreases  from  1  to  3/4  following  a  sine-squared  variation.  Furthermore,  the 
initial  field  strength  at  the  range  where  PE  computations  start  is  weighted  by  a  Gaussian 
factor  which  decays  to  -70  dB  over  the  same  altitudes.  Thus,  only  the  lowest  3/4  of  the 
FFT  results  are  actually  used  as  PE  solutions. 

To  fill  the  remainder  of  the  space,  RPO  relies  on  ray-optics  until  the  incident  ray 
to  the  earth  approaches  the  diffraction  limiL  Beyond  this  limit,  PE  results  are  utilized 
where  they  an*  avaUabie.  This  is  called  the  PE  region.  Above  the  PE  region,  if  the  space 
is  not  covered  by  ray-optics,  an  extended  optic  (XO)  region  is  defined.  The  PE  solutions 
at  the  greatest  height  are  treated  as  rays  emanating  from  below,  their  interferences  with 
the  direct  rays,  if  present,  are  taken  as  the  fields  in  this  region. 

Even  within  the  ray-optics  region,  RPO  subdivides  this  region  into  flat-earth  (FE) 
and  ray-optics  (RO)  regions.  The  earth  is  considered  as  flat  and  all  refriaction  effects  are 
ignored  in  the  FE  region,  which  is  limited  to  2.5  km  in  range  unless  the  antenna  is 
elevated  beyond  5°. 
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B.  HIGH-FREQUENCY  VERSUS  LOW-FREQUENCY  APPROXIMATION 


Bcyocid  regions  where  ray-optics  is  applicable,  RPO  relies  on  PE.  The  parabolic 
equation  approximation  to  the  Maxwell  wave  equations  is  developed  under  the  optical 
assumption  that  the  operating  frequency  is  so  high  that  a  main  direction  of  wave 
propagation  in  terms  of  a  plane  wave  can  be  specified.  Hence,  only  the  variation  of  the 
"envelope"  which  modulates  the  magnitude  and  phase  of  the  plane  wave  has  to  be 
considered,  instead  of  the  fast  variation  of  die  complete  wave  to  the  order  of  the 
wavelength.  On  the  other  hand,  M-Layer  has  its  theoretical  basis  in  an  eigenfunction 
expansion.  This  is  an  approach  most  suitable  for  low  frequency  applications  because  more 
and  more  modes  will  be  needed  as  the  size  of  the  geometrical  structure  becomes  large 
compared  to  the  wavelength.  Even  though  M-Layer  removes  the  earth  radius  as  the 
dominating  length  scale  through  the  adoption  of  a  flattened  earth,  the  thickness  of  the 
overall  atmospheric  layen  specified  is  still  far  greater  than  the  wavelength  in  all 
applications.  The  comparison  of  the  predictions  of  PE  against  M-Layer  is  one  between 
two  theories  which  approach  a  problem  from  two  extremes.  The  fact  that  the  results  of 
the  high  frequency  PE  computations  as  formulated  in  RPO  agree  almost  always  with  those 
derived  from  the  low  frequency  modal  computation  as  formulated  in  M-Layer  represents 
a  significant  engineering  achievement  by  NRaD.  However,  there  are  disagreements  at  low 
altitudes  in  the  neighborhood  of  the  radar  horizon  under  some  ducting  conditions.  These 
are  regions  of  importance  to  naval  operations.  A  defrnitive  resolution  of  these 
discrepancies  is  an  urgent  task  to  be  undertaken  in  the  immediate  future. 
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C.  SCOPE  OF  THE  WORK 


In  this  thesis,  the  radiation  of  a  vertically  polarized  transmitting  antenna  located  at 
a  height  of  IS  meters  is  investigated.  Three  modified  refi^tive  index  profiles,  with  the 
first  two  given  in  Ref.  7.  are  utilized:  a  300  m  surface-based  duct,  a  14  m  ev^oration 
duct,  and  the  combination  of  these  two  profiles  with  one  sitting  on  top  of  the  other.  Three 
frequencies  at  3  GHz,  6  GHZ  and  12  GHz  are  chosen  which  give  an  approximate 
coverage  of  the  radar  bands.  Since  only  low  altitude,  near  the  horizon  propagation  is  of 
interest,  propagation  loss  for  altitudes  up  to  100  meters  within  the  range  of  15  to  1 10  km 
are  evaluated  by  both  RPO  and  M-Layer  for  comparison. 


TL  PROGRAM  SETUP 


A.  TEST  CONSIDERATIONS 

Of  interest  to  this  thesis  research  are  the  propagation  loss  predictions  made  by  the 
parabolic  equation  approximation  method  and  by  waveguide  mode  computations.  The 
NPS  version  of  M*Layer  is  a  waveguide  mode  computation  program  which  has  been 
greatly  enhanced  in  its  accuracy  and  efficiency  [Ref.  S]  from  its  original  incarnation 
assembled  by  NOSC  [Ref.  4].  As  explained  in  Cluq)ter  I,  the  RPO  program,  completed 
by  NRaD  and  available  to  NPS  only  recently,  includes  several  methods  of  computation 
applied  to  different  regions  of  space.  To  ensure  that  only  predictions  made  under  the 
parabolic  equation  iq)proximation  are  included  in  the  comparison,  the  RPO  FORTRAN 
source  code  [Ref.  7]  is  used.  This  source  co^  consists  of  a  set  of  subroutines.  The  output 
for  each  region  comes  from  a  distinct  subroutine.  By  instructing  every  RPO  subroutine 
which  computes  the  propagation  loss  to  associate  a  code  number  to  its  output,  the 
predictions  made  with  the  parabolic  equation  approximation  are  easily  identified. 

The  RPO  FORTRAN  source  code  comes  without  a  main  program.  Nor  does  it 
contain  input  routines  and  files.  Thus,  the  first  task  of  this  research  was  to  write  a  main 
program  for  RPO  which  will  call  the  relevant  subroutines  to  compute  the  desired  ouq>ut 
Next,  data  files  had  to  be  created  to  provide  the  program  with  refiractive  index  profiles 
and  other  necessary  parameters.  An  input  sutmutine  had  to  be  written  to  feed  these  data 
to  the  main  program  and  the  RPO  sutxoutines.  Furthermore,  to  ensure  that  botii  RPO  and 


M-Laycr  use  the  same  environmental  data  and  compatible  input  parameters,  a  new  M- 
Layer  input  subroutine  was  created  which  replaces  the  subroutine  of  the  same  name  in 
the  NPS  version  of  M-Laycr.  Finally,  the  output  data  are  presented  graphically  using  a 
plotting  routine  written  in  MATLAB,  version  3.  Specific  considerations  put  into  these 
subroutines  are  discussed  in  the  following  sections,  with  the  RPO  main  program 
"RPOmain,"  the  RPO  input  routine  "RPOstdin,"  the  substituting  M-Laycr  input  routine 
"MLstdin,"  listed  respectively  in  Sections  1  through  4  of  Appendix  A. 

B.  RPO  MAIN  PROGRAM 

In  RPO,  the  parabolic  equations  are  solved  with  the  split-step  algorithm.  As  a 
prerequisite,  RPO  has  to  estimate  the  electric  field  at  an  initial  range.  It  then  proceeds, 
at  a  pre-determined  step  size,  toward  the  range  where  a  propagation  loss  prediction  is 
called  for,  computing  the  electric  field  at  every  step  along  the  way.  Hence  the  subroutines 
of  RPO  can  be  separated  into  two  types;  those  providing  other  subroutines  with  initial 
values  and  defining  their  boundaries  of  application,  and  those  handling  the  PE 
computation  at  each  step.  After  the  tedious  but  necessary  boolckeeping  of  declaring  global 
variables  and  setting  up  common  blocks  to  pass  along  values,  the  main  program  first 
reads  in  data  by  calling  the  input  subroutine  IG^tdin"  and  initializes  the  subroutines 
by  calling  "RPOini"  before  entering  into  a  loop  to  step  the  computation  through  the 
desired  range  and  write  the  propagation  loss  predictions  within  the  selected  region  to  an 
output  file.  A  flow  chart  for  the  main  program,  RPOnuun,  is  shown  in  Figure  1.  A  listing 
of  the  program  is  included  in  Appendix  A. 
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C.  INPUT  FILES  AND  SUBROUTINES 


RPO  and  M-Layer  have  different  capabilities  and  require  different  types  of  input 

For  example,  RPO  needs  a  presumed  radiation  pattern  to  set  up  its  initial  electric  field  to 

begin  computation;  it  cannot  compute  electromagnetic  (EM)  radiation  directly  from  a 

specified  current  distribution.  On  the  other  hand,  M-Layer  computes  the  radiation  from 

a  dipole  source;  it  thus  can  not  simulate  antennas  arbitrarily.  RPO  is  able  to  handle 

multiple  refractivity  profiles  at  different  ranges  due  to  further  approximations  made  in  the 

0 

theory  beyond  the  reduction  of  the  wave  equation  to  a  parabolic  equation  while  M-Layer 
is  restricted  to  deal  with  a  single  profile.  Nevertheless,  they  use  many  parameters  in 
common.  Since  the  purpose  of  this  worlc  is  to  compare  the  predictions  of  these  two 
programs,  it  is  desirable  to  have  them  read  in  exactly  the  same  files  whenever  possible. 
This  is  especially  true  for  the  refractivity  profiles. 

For  this  thesis,  three  ASCII  files  are  set  up:  one  contains  the  parameters  common 
to  both  RPO  and  M-Layer,  including  the  modified  refractivity  profile;  another  contains 
RPO  specific  paranwters  such  as  the  antenna  pattern;  the  third  contains  parameters  used 
only  by  M-Layer  stich  as  the  parameter  aloss,  which  specifies  the  greatest  range 
attenuation  rate  of  the  modes  to  be  searched,  in  dB  per  kilometer.  The  parameters  which 
are  common  to  both  RPO  and  M-Layer  are  listed  in  Table  HI.  Those  specific  to  RPO 
are  listed  in  Table  IL2.  The  M-Layer  specific  parameters  are  listed  in  Table  n.3.  Two 
batch  files  are  written  to  combine  automatically  the  two  proper  files  for  input  into  RPO 


and  M-Layer  respectively. 


Table  II.  1.  Input  parameters  common  to  RPO  and  M-Layer. 


RPO 

M-Layer 

Description 

nlevls 

nzlayr 

number  of  profile  levels^*^^ 

wind 

wind 

surface  wind  speed 

fmhz 

fqmzin 

operating  frequency 

ipolar 

mpol 

antenna  polarization^^*^^ 

ztran 

ztinit 

transmitter  height 

mlxout 

nx 

number  of  receiver  ranges 

xinit 

xinit 

inidal  range  of  the  receiver 

delx 

delx 

receiver  range  increment 

nzout 

nzr 

number  of  receiver  heights 

zrinit 

zrinit 

initial  hei^t  of  the  receiver 

delzr 

delzr 

receiver  height  increment 

zprof 

zi 

heights  at  which  profile  data 

are  specified  (an  array) 

capm 

zim 

profile  data  (modified  index  of 

refraction;  an  array) 

(a)  RPO:  from  1  to  nlevly,  M-Layer  from  0  to  nzlayr. 

(b)  follows  M-Layer  definition,  adjust  RPO  during  input 

(c)  RPO:  l=hoiizontal;  2=vertical;  M-Layen  O=horizontal;  l=vertical. 
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Table  n.2.  RPO-specific  input  parameters. 


RPO 

Description 

Remarks 

nxout 

number  of  ouq)ut  ranges  computed 

calculated  by  Xmax/delx 

where 

Xmax^XinitHlelX*mlxout 

selx 

minimum  range  to  output  data 

xprof 

ranges  at  which  index  of  refraction 

profiles  are  specified 

set  to  0.0 

nprofs 

number  of  profiles  speciHed 

set  to  1 

ipatm 

antenna  pattern 

1:  omni-direction  2:  sin(x)/x 

3:  Gaussian  4:  cosecant-squared 

5:  height-finder  6:  user  defined 

use  1  or  3  only 

beamw 

antenna  elevation  beamwidth 

elang 

antenna  elevation  angle 

set  to  0.0 

rfacs 

number  of  height-finder  data 

set  to  0  (not  used) 

iscatt 

includes  troposcatter 

0:  no  1:  yes 

set  to  0  (not  used) 

Table  n.3,  M-Laycr-spccific  parameters. 


M-Layer 

Description 

Remarks 

mfile 

0:  read  input  and  compute  eigenvalues 

1:  read  eigenvalues  as  input 

delfq 

frequency  increment 

set  to  0.0  (not  used) 

r^eq 

number  of  frequencies  to  be  used 

set  to  1 

aloss 

maximum  range  attenuation  rate  in 

set  to  2.0  or  S.O 

dB/km  of  modes  to  be  found 

delzt 

height  increment  of  transmitter 

set  to  0.0  (not  used) 

nzt 

number  of  transmitter  heights 

set  to  1 

r«fi 

reference  height  at  which  refht  and 

set  to  0.0 

rifgab  are  given 

refm 

modified  refractivity  at  the 

set  to  339.0 

reference  height 

rtfgab 

modified  gas  absorption  at  the 

set  to  0.0  (not  used) 

reference  height 
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In  RPO,  special  formulas  for  the  dielectric  constant  and  conductivity  of  sea  water 
as  a  function  of  frequency  are  used.  Sea  surface  roughness  is  also  given  in  terms  of  a 
function  of  wind  speed.  To  incorporate  this  into  M-Layer,  the  input  subroutine  "MLstdin" 
for  M-Layer  is  modified  to  cany  out  these  computations  using  the  same  formulas  before 
providing  these  parameters  to  the  program,  even  though  surface  roughness  is  not 
considered  in  this  comparison. 

D.  VALIDATION 

Before  proceeding  with  the  comparison,  a  few  RPO  test  cases  are  run  using  this 
newly  written  main  program  and  the  input  subroutine.  Specifically,  the  tests  listed  in 
Ref.  7  under  the  names  LOBW  (low  beam  width  limit),  SBDUCT  (surface-based  duct) 
and  EDUCT  (evaporation  duct)  are  carried  out  The  tabulated  results  in  Ref.  7  are 
reproduced  exactly  as  long  as  the  range  increment  detx  and  height  increment  delz  are  the 
same  as  those  specified  therein.  Otherwise,  variations  of  up  to  4%  are  observed. 


m.  RESULTS 


Propagation  of  waves  through  several  refractivity  profiles  at  many  frequencies  has 
been  investigated.  A  clear  trend  has  emerged  which  shows  disagreement  between  the 
results  from  RPO  and  M-Layer  at  ranges  near  the  horizon  when  surface  based  ducts  are 
involved.  In  this  thesis,  results  from  three  profiles  are  presented:  a  300  m  surface-based 
duct  which  is  specified  and  used  in  Ref.  7  under  the  test  name  SBDUCT;  a  14  meter 
evaporation  duct  which  is  also  specifred  and  used  in  Ref.  7,  under  the  test  name  EDUCT; 
and  a  combination  of  these  two  ducts  by  merging  the  profile  of  the  evaporation  duct  with 
that  of  the  surface-based  duct  The  modified  refractivity  profile  of  the  300  m  surface- 
based  duct  is  given  in  Table  in.l.  The  profile  for  the  14  m  evaporation  duct  is  given  in 

Table  nU.n«proaie  for  the  combin«ion  Table  in.l.  A  300  m  sorfece-based  duct 
is  given  in  Table  in.3.  For  each  of  the 
three  profiles,  three  frequencies  at  3  GHz, 

6  GHz  and  12  GHz  are  selected.  With  the 
transmitter  fixed  at  a  height  of  IS  m,  the 
propagation  loss  of  up  to  100  m  at  ranges 
of  15,  20,  30,  40,  50,  60,  70,  80,  90,  100 
and  110  km  is  plotted.  Since  the  results  at  different  frequencies  show  similar  features, 
only  those  at  6  GHz  are  included  in  this  chapter.  Those  at  3  GHz  and  12  GHz  ate 
collected  in  Appendices  B  through  D. 


i 

Zi 

metoa 

Mi 

0 

0.000 

339.0 

1 

250.0 

368.5 

2 

300.0 

319.0 

3 

1000.0 

401.6 
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Table  III.2.  A  14  m  evaporation  duct 


Zi 

meters 

Mi 

0.000 

339.00 

0.040 

335.10 

0.100 

333.66 

0.200 

332.60 

0.398 

331.54 

0.794 

330.51 

1.585 

329.53 

3.162 

328.65 

6.310 

327.96 

12.589 

327.68 

14.000 

327.67 

25.119 

328.13 

39.811 

329.25 

50.119 

330.18 

63.096 

331.44 

79.433 

333.12 

100.000 

335.33 

125.893 

338.20 

158.489 

341.92 

199.526 

346.69 

209.526 

347.87 

Table  in.3.  A  300  m  surface-based 
duct  over  a  14  m  evaporation  duct 


Zi  Mi 

meters 


0.000  339.00 


0.040  335.10 


0.100  333.66 


0.200  332.60 


0.398  331.54 


0.794  330.51 


1.585  329.53 


3.162  328.65 


6.310  327.96 


12.589  327.68 


14.000  327.67 


25.119  328.13 


39.811  329.25 


50.119  330.18 


63.096  331.44 


79.433  333.12 


100.000  335.33 


125.893  338.20 


158.489  341.92 


199.526  346.69 


209.526  347.87 


250.0 


300.0  319.0 


1000.0  401.6 


For  all  the  cases,  the  polarization  is  chosen  to  be  vertical.  The  receiver  height 
increment  delzr  is  set  at  0.5  m.  The  receiver  range  increment  delx  is  set  at  2500  m.  For 
RPO  computations,  the  maximum  height  zmax  is  set  at  100  m;  the  maximum  range  xmax 
is  set  at  115  km.  The  Gaussian  beam  pattern  with  a  beamwidth  of  5°  is  chosen  for  the 
antenna.  Within  the  altitudes  and  ranges  considered,  there  is  no  perceivable  difference 
when  the  omni-directional  pattern  is  used  instead.  For  M-Layer  computation,  the 
parameter  ahss  is  set  to  2  dB/km.  In  all  the  figures,  results  from  M-Layer  are  drawn  as 
a  solid  line.  Results  from  RPO  are  marked  with  asterisks  if  they  are  in  the  PE  region  and 
with  dots  if  they  are  in  the  RO  region.  For  easy  reference,  the  height  of  the  transmitter 
horizon  at  each  distance,  based  on  the  four-thirds  effective  earth  radius,  is  indicated  with 
a  horizontal  line  drawn  across  the  figure  on  which  it  is  present  It  starts  to  ^pear  in  the 
figure  for  the  20  km  range.  Figures  beyond  50  km  lie  below  the  horizon  completely  and 
this  horizontal  line  cannot  be  seen. 

In  what  follows,  results  of  RPO  and  M-Layer  computations  at  6  GHz  are  presented. 
Their  analysis  and  discussions  are  given  in  Chapter  IV. 

A.  300  METER  SURFACE-BASED  DUCT 

Figures  2  through  12  show  the  propagation  loss  at  various  ranges  from  15  through 
1 10  km  when  a  300  m  surface-based  duct  is  present  This  is  the  profile  in  which  most 
signifreant  deviations  between  RPO  and  M-Layer  are  observed,  especially  at  50  and  60 
km  ranges.  Even  at  40  km,  the  two  differ  by  10  to  30  dB  at  low  altitudes.  On  the  other 
hand,  the  two  agree  well  within  line-of-sight  and  deep  shadow  regions. 
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Figure  3.  Propagation  loss  at  20  kia 
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Figure  5.  Propagation  loss  at  40  km. 
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Figure  8.  Propagation  loss  at  70  km. 


Figure  9.  Propagation  loss  at  80  km. 
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Figure  10.  Propagation  loss  at  90  km. 


Figure  11.  Propagation  loss  at  100  km. 
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Figure  12.  Propagation  loss  at  1 10  km. 


B.  14  METER  EVAPORATION  DUCT 

Figures  13  through  23  show  the  propagation  loss  in  the  presence  of  a  14  m 
evaporation  duct  RPO  and  M-Layer  agree  well  over  the  entire  range.  From  the  figures 
in  Appendix  C  ,  it  can  be  seen  that  this  is  true  for  all  the  frequencies  investigated,  with 
only  a  less  than  1.5  dB  difference  around  where  the  loss  is  maximum  for  the  12  GHz  case 
in  the  over-the-horizon  region. 
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Figure  13.  Propagation  loss  at  15  km. 
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Figure  15.  Propagation  loss  at  30  km. 


Figure  16.  Propagation  loss  at  40  km. 
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Figure  18.  Propagation  loss  at  60  km. 
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21.  Propagation  loss  at  90  km. 


Figure  22.  Propagation  loss  at  100  km. 
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Figure  23.  Propagation  loss  at  1 10  km. 

C.  SURFACE-BASED  DUCT  OVER  EVAPORATION  DUCT 

The  agreement  in  one  case  and  the  disagreement  in  another  between  RPO  and  M- 
Layer  prompted  an  investigation  into  a  combined  profile  consisting  of  these  two  ducts. 
Figures  24  through  34  show  Ac  propagation  loss  under  such  a  refiractivity  profile.  At 
6  GHz,  RTO  and  M-Layer  agree  well  over  Ac  entire  range  except  between  50  and  60  km 
when  M-Layer  displays  an  mcrcasc  m  field  strengA  Award  Ac  surface  over  Ac  lowest 
few  meters  while  RPO  continues  to  decrease.  This  leads  to  more  Aan  10  dB  deviations 
in  propagation  loss  at  Aesc  ranges.  At  12  GHz,  Ac  two  programs  agree  even  better  when 
M-Layer  docs  not  show  Ac  mcrease  toward  surface  level  within  Ac  last  few  meters.  At 
3  GHz,  this  diffCTcnce  is  more  pronounced.  Below  10  to  20  m,  Ais  Avcrgence  starts  to 
appear  and  leads  to  deviations  in  propagation  loss  of  up  to  40  dB. 
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Igure  25.  Propagation  loss  at  20  km. 
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i  26.  Propagadon  loss  at  30  km. 


Figure  27.  Propagation  loss  at  40  km. 
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t  28.  Propagation  loss  at  50  km. 


Figure  29.  Propagation  loss  at  60  km. 
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e  30.  Propagation  loss  at  70  km. 


Figure  31.  Propagation  loss  at  80  km. 
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Figure  33.  Propagation  loss  at  100  km. 
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IV.  ANALYSIS  AND  CONCLUSIONS 


A.  ANALYSIS 

Since  RPO  is  a  high  frequency  approximation  while  M-Layer  is  based  on  a  low 
frequency  technique,  and  the  vertical  extent  within  which  the  index  of  refraction  of  the 
atmosphere  is  specified  is  much  larger  than  the  wavelength  under  consideration,  RPO  is 
expected  to  be  more  accurate  at  short  ranges.  The  fact  that  M-Layer  matches  the  results 
of  ray-optics  computations  in  Fig.  13  at  15  km  from  the  transmitter  for  heights  over  90 
m  above  ground  confirms  the  reliability  of  M-Layer.  RPO  and  M-Layer  agree  well  also 
in  the  region  far  beyond  the  radar  horizon.  What  happens  then  over  the  region  in-between 
when  a  300  m  surface-based  duct  is  present? 

The  physics  of  the  situation  is  revealed  if  the  earth-flattening  approximation  as 
prescribed  in  Ref.  6  is  re-examined.  In  fact,  both  RPO  and  M-Layer  are  treating  tiie 
propagation  problem  as  one  involving  only  a  flat  "earth"  surface.  The  waves  are  traveling 
in  a  flat,  layered  dielectric  waveguide.  There  is  no  blocka^  of  rays  due  to  the  earth  to 
create  a  shadow;  instead,  the  rays  are  bent  away  from  the  "shadow"  region  by  the  linearly 
increasing  component  of  the  modified  refractivity.  The  ducting  structure  bends  the  wave 
back  towards  the  surface.  Properly  implemented,  these  two  programs  should  provide 
identical  predictions  over  ranges  of  common  validity. 

Above  a  flat  surface,  the  concept  of  interference  between  a  direct  ray  and  a  reflected 
ray,  even  though  bent  by  the  atmosphere,  remains  valid.  The  bending  shifts  the  locations 
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of  constructive  and  destructive  interferences,  but  their  separations  arc  less  affected.  From 
Figs.  2  and.  13,  it  is  clear  that  at  IS  km,  the  separations  between  two  neighboring  nulls 
above  ground  are  about  28  m,  independent  of  the  particular  environment  and  close  to  a 
free-space,  parallel-ray  estimate  of  2S  m.  The  separations  of  nulls  are  expected  to  increase 
as  the  range  is  increased.  This  is  observed  in  Rgs.  13  through  23  for  the  propagation  in 
the  presence  of  the  14  m  evaporation  duct  But  when  the  300  m  surface-based  duct  is 
present  both  RPO  and  M-Layer  predict  for  ranges  greater  than  90  km,  similar  field 
strength  variations  which  oscillate  much  faster  than  those  observed  within  the  line-of-  , 
sight  region.  It  is  clear  that  this  is  a  phenomenon  due  to  the  presence  of  the  surface-based 
duct  In  terms  of  high  frequency  ray-optics  terminology,  this  fast  oscillation  is  due  to 
additional  rays  which  are  bent  back  by  the  duct  In  terms  of  waveguide  mode  theory,  it 
is  due  to  the  propagation  of  waves  in  many  modes  established  in  the  surface-based  duct, 
each  at  a  distinct  phase  velocity,  thus  interfering  severely. 

The  devi^ons  between  RPO  and  M-Layer  at  ranges  from  30  through  70  km  which 
show  up  in  Figs.  4  to  12  can  now  be  explained.  RPO  fails  to  fully  take  into  account  the 
effects  of  a  high  surface-based  duct  This  argument  is  further  supported  by  investigating 
the  effects  of  lowering  the  duct  height  At  all  frequencies,  the  fast  oscillation 
characteristic  of  the  presence  of  a  surface  based  duct  sets  in  at  a  shorter  range  with  a 
lower  duct  height  Agreement  between  RPO  and  M-Layer  gets  better  as  the  duct  height 
is  lowered.  These  can  be  observed  from  the  results  presented  in  Appendix  E  for 
propagation  in  a  150  m  surface-based  duct  and  in  Appendix  F  for  propagation  in  a  100  m 
duct 
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B.  CONCLUSIONS 

At  low  altitudes  in  the  neighborhood  of  the  radar  horizon,  deviations  between  RPO 
and  M-Layer  can  be  significant  RPO  appears  not  to  be  able  to  properly  include  the 
effects  of  a  high  surface-based  duct  at  a  short  range.  It  is  recommended  that  parabolic 
equation  computations  should  start  at  a  closer  range  than  currendy  prescribed  in  RPO. 
The  altitude  at  which  the  filtering  of  field  strength  starts  should  always  be  much  higher 
than  the  duct  height 
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APPENDIX  A;  PROGRAM  SOURCE  CODES 
This  Appendix  contains  listings  of  the  RPO  main  program  RPOmain,  the  RPO  input 
subroutine  RPOstdin,  the  M-Layer  input  sutxoutine  MLstdin  and  the  MATLAB  M-file 
RPOMLA.M  for  plotting  the  propagation  loss  computed  with  RPO  and  M-Layer. 

1.  Subroutine  RPOmain 


program  RPOmain 
c 

c  PURPOSE:  Main  progr2un  for  RPO. 

c 

c 

c  input : 

c  Argvunent  List:  None 

c  Common:  fmhz,  ztran.  ipolar,  ipatrn,  beamw,  elang,  hfang,  hffac,  nfacs 
c  capm,  zprof,  nprofs,  nlevls,  xprof,  wind 

c  iscatt,  nuoclev,  maxnx,  nxout,  meocnz,  nzout,  maxpro,  xnuuc,  zmax 

c 

c 

c  output: 

c  Argument  list:  None 
c  Common:  losscb,  smg 
c 
c 

c  subroutines  calling  RPOmain:  None. 

c  subroutines  called  by  RPOmain:  rpostdin,  rpoini,  rpostp, 

c 

c 

c  common  block 


c 

/system/ 

c 

fmhz 

EM  system  frequency 

c 

ztran  : 

antenna  height 

c 

ipolar: 

euitenna  polarization 

c 

ipatm: 

antenna  pattern 

c 

beamw  : 

euitenna  vertical  beeun  width 

c 

elang  : 

emtenna  elevation  angle 

c 

hfemg  : 

height -finder  angles  array  in  degrees  (0.  to  99) 

c 

hffac  : 

height-finder  power  reduction  factor  array  (0.  to  1.0  ) 

c 

nfacs  : 

number  of  power  reduction  angles /factors  for  user-defined 

c 

height -finder  radar 

c 

c 

/enviro/ 

c 

capn 

(•)  profile  modified  refractivity  array 

c 

zprof  : 

(*)  profile  heights  in  meters  (.GE.  0.) 

c 

nprof 3 : 

number  of  profile  levels 

c 

nlevls: 

number  of  refractivity  profile  levels  (1  to  mzuclev) 

c 

xprof  : 

(*)  range  to  each  profile  in  meter  (.GE.  0.) 

c 

c 

wind  : 

wind  speed  at  range  zero  in  meter/ sec  (.GE.  0.) 
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c 

/init/ 

c 

iscatt 

c 

maxlev 

c 

maxnx 

c 

nxout 

c 

maxnz 

c 

nzout 

c 

maxpro 

c 

xmAX 

c 

zmax 

c  Declares  the 


flag  to  include  troposcatter 

meucimuin  number  of  profile  level  above  zero  height 
maximum  number  of  output  range  points 
number  of  output  range  points  (1  to  maxnx) 
meucimum  number  of  output  height  points 
number  of  output  height  points  (1  to  maxnz) 
maximun  number  of  profiles  beyond  zero  range 
maximum  range  for  output  in  meter 
nAximum  height  for  output  in  meter 
types  of  parameters 


character*8  filein 


integer*4  ipolar,  ipatrn.  nfacs,  nprofs,  nlevls,  iscatt, 
*  nxQut,  nzout.  inaxlev,  iMUcnx,  maxnz,  maxpro 

integer*2  losscb,  srng 

real  fmhz,  ztran,  beamw,  elang,  hfang,  hffac, 

>  xprof,  wind,  capm,  zprof,  xmax,  zimuc 

c  The  following  include  file  contains  a  PARAMETER  statement  to 
c  define  meocimun  array  sizes  meucpro,  maxlev,  meucnz 

INCLUDE  'RPOSIZE.INC' 


c 

c  The  following  PARAMETER  statement  defines  meucimum  array 
c  dimensions  used  throughout  RPO.  Generally,  these  constants 

c  will  have  to  be  changed  for  each  inqplementation  of  RPO. 
c  GLOSSARY ; 

c  maxlev:  nuucimum  number  of  profile  levels  above  zero  height, 

c  maxnx  :  maximum  number  of  output  range  points, 

c  maxnz  :  maximum  number  of  output  height  points, 

c  maxpro:  meucimum  nuunber  of  profiles  beyond  zero  range, 

c 

PARAMETER  (nKUcnx  a  440,  meucnz  a  280,  maxlev  a  50,  metxpro  a32} 
c 

diroens  ion  losscb  (meuoiz ) ,  srng  (nuucnz ) 

COMMON  /system/  fmhz,  ztreui,  ipolar,  ipatm,  beeunw,  eleuig, 
hfang(lO),  hffac(lO),  nfacs 

COMMON  /enviro/  nlevls,  nprofs,  wind,  zprof (0 :m2uclev,  0:naxpro), 
capm(0:maxlev,  0:nuucpro),  xprof  (0  :maxpro) 

COMMON  /init/  iscatt,  nxout,  nzout,  xmax,  zimuc 

COMMON  /misc/  jminPE,  nx,  xstep,  zstep,  wl,  rk,  fterm,  pi 

COMMON  /inout/  filein,  selx,  mlxout 


c  read  data  from  input  files 

writs  (*,*)  'Begin  with  rpostdin' 
call  rpostdin 

c  initialize  RPO 

write  (*,*)  'Begin  with  rpoini' 
call  rpoini (nsteps) 


c  call  rpostp  to  ccsnpute  propagation  loss 
c  and  write  the  result  to  a  file  on  disk 
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op«n(16, filaxfilein// ' .out') 
do  ial,  nxout 

call  rpostp(x, losscb. srng) 

.Xkmxx* .  001 

write(*,*)  'Begin  with  xpostp ' , i , xkn 
if  (x.ge.selx)  then 
do  jsl.nzout 

write  (16,1102)  xlon,  j*zstep,  .l*losscb( j) ,  srng(j) 
end  do 
end  if 
end  do 

1102  forTnat(f9.2,  6x.  f6.1.  5x.  f6.2.  5x,  i2) 
end 


2.  Subroutine  RPOstdin 


c** ****•••***•********•*******•****•***•*•*•**••*•*••**••••* •*••****•’ 

c  Subroutine  RPOstdin  reads  in  data*  from  a  file  concocted  from  two 
c  files:  one  contains  parameters  specific  to  RPO,  the  other  contains 
c  and  parameters  used  by  both  RPO  and  M- Layer,  including  the  modified 
c  refractivity,  the  treuismitter  height,  and  the  output  points. 

(;*«***•*••«***•*«•«••«•••**«*«•••*«»*»••*••*•*****•*******•••**••••••< 

subroutine  RPOstdin 

c  The  following  include  file  contains  a  PARAMETER  statement  to 
c  define  m2ueimun  array  sizes  maxpro,  meuelev 
INCLUDE  'RPOSIZE.INC' 
c 

c  The  following  PARAMETER  statement  defines  maximum  array 
c  dimensions  used  throughout  RPO.  Generally,  these  constants 

c  will  have  to  be  changed  for  each  in^lementation  of  RPO. 
c  GLOSSARY : 

c  maxlev:  meucimuro  number  of  profile  levels  eUbove  zero  height, 

c  maxnx  :  maximum  number  of  output  range  points, 

c  maxnz  :  maximum  number  of  output  height  points, 

c  maxpro:  maximum  number  of  profiles  beyond  zero  range, 

c 

PARAMETER  (maxnx  s  440,  maxnz  s  280,  maxlev  s  50,  meucpro  s32) 
c 

character*8  filein,  fileb 

COMMON  /system/  fmhz,  ztran,  ipolar,  ipatrn,  beeunw,  eleuig, 

*  hfeuig(lO),  hffac(lO),  nfacs 

COMMON  /enviro/  nlevls,  nprofs,  wind,  zprof (0 :maxlev,  Ormaxpro), 
capm(0:maxlev,  0:maxpro),  xprof (0:inaxpro) 

COMMON  /init/  iscatt,  nxout,  nzout,  xmax,  znux 
COMMON  /inout/  filein,  selx,  mlxout 

c - read  the  RPO  parameter - 

read(*, ' (a) ' )  filein 
read ( * , * )  nxout 
read(*,*)  selx 
read ( * , * )  xprof ( 0 ) 
read ( * , * )  nprofs 

c  RPO  starts  its  array  index  with  0.  It  adopts  Microsoft  BASIC 
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c  convention  in  its  coding  of  this  FORTRAN  prograun.  In  the  documentation 
c  and  the  input  file,  'nprofs*  is  the  total  number  of  profiles.  Thus 
c  *nprofs*  as  the  array  index  as  used  in  RPO  has  to  be  adjusted  by 
c  substracting  "nprofs*  by  1. 
nprofs«nprofs>l 
read(*,*)  ipatrn 
read  ( * ,  * )  beeunw 

c  For  omnidirectional  pattern,  the  parameter  *be«uiiw*  is  ignored  in  the 
c  progreun.  On  the  other  hand,  the  upper  limit  for  'beamw*  is  45  degrees, 
if  ( ipatrn. eq.l)  then 
beainwx45.0 
end  if 

read(*,*)  elang 
read(*,*)  nfacs 
read{*,*)  iscatt 

c - read  parameters  common  to  RPO  ^md  M-Layer - 

read(*, ' (a) ’ )  fileb 
read(*,*)  nlevls 

read { * , * )  wind  , 

read ( * , * )  fmhz 
read(*,*)  ipolar 

c  The  definition  of  ipolar  in  RPO  is  1  plus  that  in  M~LAYER.  ' 
ipolarxipolar-t-l 
read(*,*)  ztrem 
read(*,*)  nzout 
raad{*,*)  delzr 
read(*,*)  zrinit 
zmeuc3nzout*delzr>zrinit 
read ( * , * )  mlxout 
read{*,*)  delx 
read ( * , • )  xinit 
delxsdelx^lOOO 
xinitaxinit*1000 
xnmxxmlxout  *de lx>xini t 
do  ixO, nlevls 

read ( * , * )  zprof ( i , 0 ) 
read ( * , * )  capm ( i , 0 ) 
capm(i,  0)  3capm(i,  0)  -1-339 .0 
end  do 

open(17, filexfilein// ' . in' ) 

write(17,*)' - parameters  coinnon  to  RPO  and  M-Layer - ' 

write (17, *) filein,  '  :RPO  file' 

wr ite ( 17, *) nlevls,  '  :nlevls  (count  from  0)' 
write (17, 190) wind 
190  format(fi5.2,  '  :wind') 
write ( 17 , 192 ) fmhz 
192  format(fl5.2,  '  ifmhz') 

wr ite ( 17, *) ipolar- 1,  '  :ipolar  (0:horizontal  Irvertical  )' 
write ( 17 , 194 ) ztran 
194  format (f 15 .2,  '  : ztran') 

write(17, ♦) nzout,  '  tnzout' 

write (17,201) delzr 

201  format(fl5.2,  '  : delzr') 

write (17,202) zrinit 

202  format (f 15 .2,  '  : zrinit’) 

write(17, *) mlxout  ,'  :mlxout  (nuuc  number  of  remge  output)  ' 

write (17,203) delx 

203  format (f 15. 2,  '  :delx') 


write (17, 204)xinit 

204  format ( CIS .2,  '  ;xinit  *) 

write{17,*)' -  RPO  parameters  - ' 

write (17, 205) zmeuc 

205  fonnat(fi5.2,  '  :zmax  (nzout*delzr, calculated  by  program)'  ) 
write ( 17 , 206 ) xmax 

206  format(fl5.2,  '  :xmax  (xinit4-delx*mlxout, calculated  by  program#)  ' ) 

write(  17 ,  *) nxout ,  '  :nxout  (xinauc  /  delx,  read  from  input)  ' 

write ( 17 , 207 ) selx 

207  format ( CIS .2,  '  :selx  (selected  range  x  >s  selx  )') 
write (17, 208)xproC(0) 

208  Cormat(fl5.2,  '  ixprof') 

write ( 17, •)nprofs+l,  '  snprofs' 
write ( 17, •) ipatrn,  '  tipatrn' 

write  (17, 209)beaunw 

209  format(£15.2,  '  :beaunw') 

write (17,210) elang 

210  format ( f 15 .2,  '  :elang') 

write(17,*)  '  N/A  :h£ang  (not  used)' 

write(17,*)  '  N/A  :h££ac  (not  used}' 

write ( 17 , *)n£acs,  '  :n£acs' 

write(17, •) iscatt,  '  tiscatt' 

write (17,*)' - M-pro£ile - value  on  sur£ace:  339.0  - ' 

do  i=0,nlevls 

write{17, 211) zpro£{i, 0) 

211  £ormat (£15 .2,  '  :zpro£') 
write (17,212) capm  ( i , 0 ) 

212  £ormat(£15.2,  '  :capm  ') 
end  do 

return 

end 


3.  Subroutine  MLstdin 


^««****«*****«****«**«««****«****««*««**«****«*«**««*««*«««****«««***«« 
c  Subroutine  MLstdin  reads  data  £rom  a  £ile  concocted  £rom  two  £iles:  * 
c  one  contains  parauneters  speci£ic  to  M-Layer,  the  other  contains  * 

c  parauneters  used  by  both  RPO  and  M-lrfiyer,  including  the  modi£ied  * 

c  re£ractivity,  the  transmitter  height,  and  the  output  points.  * 

(.*********•**«**•«*•«*•«*«««*«««*•**•«***•******•••**«*****«•••••**•«** 

subroutine  MLstdin 
c 

c  MLstdin  is  the  revised  input  prograus  o£  M-Layer,  NPS  version,  to 

c  read  in  data  £iles.  The  common  block  /inpt9/  has  been  removed  £rom 

c  all  subroutines, 
c 

implicit  real*8  (a-h,o-z) 
complex* 16  qeigen 

integer  i£lgab,mpol,nzt,nzr,nx,nzlayr,nrroode,m£ile, i 
character*40  £ilein,  £ileb 

c 

c  use  include  £ile  £or  parauneters  o£ 
c  mxlayr  mauc  #  layers 

c  mxmode  roauc  #  modes 

c 

include  ' mlaparm . inc ' 


42 


c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


c 


include  £ile  to  define  the 

maximum  #  of  layers  (mxlayr) 
maximum  #  of  modes  (mxmode) 


par2uneter  (mxlayrs35  ) 
pareuoeter  (mxmodes390) 

dimension  zi  (mxlayr -t-l) .  zim(mxlayr-^l) ,  zigab{mxlayr-»'l) , 
$  qeigen (mxmode) 


COTinon 

+ 

♦ 

+ 


/inptO/filein,m£ile,  fcpnzin.nqpol.aloss.dielcg,  sigmag, 

ztinit . delzt , nzt. zrinit. delzr, nzr, xinit , delx, nx, 
ref  z,  refm,  refgeds,  zim,  zigab 
/inptl/nzlayr 
/inpt2/zi 
/inpt4/rmsbht 
/modes/nrmode, qeigen 


read  M-LAYER  specific  parameters 
read( ' (a) ' )  filein 
raad(*,*)  mfile 
read(*,*)  delfq 
read(*,*)  nfreq 
read(*,*)  aloss 
read(*,*)  iflged> 
read(*,*)  delzt 
read { * , * )  nzt 
read(*,*)  zref 
read(*,*)  refz 
read(*,*)  refm 
read(*,*)  refgab 

read  parameters  common  to  RPO  and  M^Layer 
read( ' (a) ' )  fileb 
read(*,*)  nzlayr 
read{*,*)  wind 
read(*,*)  fqmzin 
read(*,*)  mpol 
road{*,*)  ztinit 
read ( * , * )  nzr 
read(*,*)  delzr 
read(*,*)  zrinit 
read ( * , * )  nx 
read(*,*)  delx 
readj*,*)  xinit 


c  The  profile  must  contain  at  least  three  levels.  The  M  gradients 

c  in  adjacent  layers  must  not  be  equal, 

c 

do  isl,nzlayr+l 
read{*,*)  zi(i) 
read(*,*)  zim(i) 
end  do 

if  (mfile. ne.O)  then 

c  mfileoO  indicates  input  file  contains  eigenvalues 
read{*, *)nrmode 
do  isl,nrmode 
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read(*,*)  qeigen(i) 

•nd  do 
•ndif 

c  Calculate  the  root  mean  square  bump  height  of  sea  surface. 

c 

nnsbhtsO . 00514*wind**2 
c 

c  Calculate  the  dielectric  constant  and  conductivity  of  sea  water. 

if  (fcipnzin  .LE.  1500.)  then 
dielcgs  80. 
sigma  3  4.3 

else  if  (fqmzin  .LE.  3000.)  then 

dielcgs  80.  -  .00733  •  (fqmzin-1500 . ) 
sigma  3  4.3  .00148  *  (fqmzin-1500.) 

else  if  (fqmzin  .LE.  10000.)  then 

dielcg=  69.  -  .00243  *  ( fqmzin-3000 . ) 
sigma  «  6.52  +  .001314  *  {f<jnzin-3000 . ) 
else 

dialcg3  51.99 
sigma  3  15.718 
end  if 

c - write  input  data  to  a  file  on  the  disk - 

open(27, f ilesf ilein// ' .in' ) 

write (27,*)  ' -  parameters  common  to  RPO  and  M-Layer  - 

write(27,'  (2(A13) ) ' ) filein,  ' :ML  file' 
write ( 27, *)nzlayr,  '  :nzlayr' 

write (27, 190)wind 
190  format (f 15 .2,  '  :wind') 
write (27 , 192) fqmzin 
192  format(fl5.2,  '  :fqmzin') 

write(27,  *)mpol,  '  :inpol  (0:horizontal  l:vertical  )' 

write{27 , 194) ztinit 
194  format(fl5.2,  '  :ztinit') 

write ( 27, *)nzr,  '  ;nzr' 

write(27,201)delzr 

201  format(fl5.2,  '  :delzr') 
write (27,202) zrinit 

202  format(fl5.2,  '  : zrinit') 

write (27, *)nx,  '  ;nx' 

write (27,203) delx 

203  format (f 15. 2,  '  :delx  in  km') 
write (27, 204) xinit 

204  format (f 15. 2,  '  :xinit  in  km') 

write(27,*)  ' -  M-Layer  specific  parameters  - 

write{27,*)  mfile,  '  :mfilesO:read  input  and  compute  e.v.' 

write (27, 205 )delfq 

205  format (f 15. 2,  '  rdelfq') 

write(27, *)nfreq, '  ;nfreq' 

write (27,207) aloss 

207  format (f 15. 2,  '  raloss') 

write(27, *) iflgab, '  riflgab' 

write (27, 208) delzt 

208  format (f 15. 2,  '  ;delzt') 

write(27, *)nzt,  '  :nzt' 
write(27, 209) zref 

209  format(fl5.2,  '  :zre£') 
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writ«(27, 210) r«£z 


210 

format (f  15. 2.  ' 

:refz‘ ) 

write(27, 211)refm 

211 

format (f 15. 2,  ' 

:refm' ) 

write(27. 212)  refgab 

212 

forraat(fl5.2,  ' 

:  refgab' ) 

write(27, *) ' - 

—  M-profila  - 

do  ial.nzlayr-fl 

%irrite(27,214) 

zi(i) , i 

214 

format (f 15. 2, 

'  :2i('  .12. 

write (27, 215) 

zim(i) . i 

215 

format (f 15. 2, 

'  :zim('  .12, 

end  do 

return 

end 

4.  MATLAB  Plotting  M-file 


%  file  name: *RPOMLA.M 
clg 

clear  ' 

£nosinput('  input  the  no.  of  test  file: 

%  delete  the  old  .met  file 

xdels( 'delete  rl'.  fno,  '.met']; 
eval (xdel) ; 

%  load  the  data 

tls[ 'load  d:\matlab\wu\dat\' , 'rpo',  fno,  '.dat'l 
t23 [ ' load  d: \matlab\wu\dat\ ' , 'ml' ,  fno,  '.dat']; 
eval (tl) ; 
eval(t2) ; 

%  read  the  data 
rpfr»['rpo',  fnoj; 
mlfrsi'ml',  fno] ; 

rl=Irpfr, r2slrpfr, ' { : ,2) ' ] ; 

r3a[rpfr, ' { :,3) ' ]  ;  r4=[rpfr, '(:,4)']; 
mlsCmlfr, ' ( ; , 1) ' ] ;  m2a(mlfr. ' ( : , 3) ' j ;  m3s[mlfr, ' ( : , 6) ' ] ; 
rpoxaeval ( r 1 ) ;  rpozaeval ( r2 ) ;  rposseval ( r3 ) ;  rporseval ( r4 ) ; 
mlxseval (ml ) ;  mlzseval (m2 ) ;  mlsseval  (m3 ) ; 

%  set  up  initial  conditoin  for  «diile  loop 

euiacl; 

profsinput(  "The  name  for  profile:  ','s'); 
frqsinput ( "The  frequency  (in  MHz):  '); 
anhs input ( 'The  antenna  height  (m) :  '); 

%  loop  for  range  increment 
while  em>0 

xrg«  input ,  '  The  x  range  '  ) ; 

%  calculate  the  receiver  antenna  height  of  horizontal  disteince 
aes€371000; 

eunh2a  (  xrg*1000-sqrt  (2*ae*anh)  )''2  /(2*ae); 

)cl«l;  k2>l;  Ic3sl;  k4sl;  p«l; 

%  xl,  x2...  is  used  to  decide  which  region  should  be  plotted  in  plot. 
xl«0;  x230;  xSsO;  x430; 

%  find  the  index  of  rpo  loss  for  x  range  equal  to  xrg 
tmpxrg«xrg-«- .  05 ; 

xindaxsfind(  (rpox>axrg)  &  (rpox<atnpxrg)  ) ; 
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3tri=xindex{ 1) ;  endi=xindex(  length (xindex)  ) ; 

%  find  the  min.  <ind  max.  of  rpo  loss  for  x  range  equal  to  xrg 
xrpgjninsminCrpos  (xindex) )  ; 
xrpom^LX3max(^pos  (xindex)  )  ; 

%  get  rpo  loss,  receiver  height,  according  to  the  method  of  computation 
%  (fe,  ro,  xo  &  pe) 

for  isstritendi 

if  rpor(i)*=l 

rpzl()cl)=  rpoz(i);  rpsl(lcl)=  rpos(i);  rprl()cl)=  rpor(i); 

)tl=Jcl-t-l; 

xl=l; 

end 

if  ’•por(i)==2 

rpz2(k2)3  rpoz(i);  rps2(lc2)s  rpos(i); 

k2=)c2+l; 

x232; 

end 

if  rpor(i)==3 

rpz3(lc3)=  rpoz(i);  rps3(k3)=  rpos(i); 

Ic3=:)t3  +  1; 

X3=3;  i 

end 

if  rpor(i)==4 

rp24{k4)=  rpoz(i);  rps4(k4)=  rpos(i); 
k43k4-^l; 

X4s4; 

end 

end 

%  find  the  index  of  mlayer  loss  for  the  x  range  equal  to  xrg 
xindexsfind(  (mlx>sxrg}  &  (mlx<stn^xrg)  } ; 
strisxindexU) ;  endi3Xindex(  length(xindex)  ); 

%  find  the  min.  and  nuuc.  of  ml  loss  for  x  range  is  equal  to  xrg 
xmlminsmin (mis (xindex) ) ; 
xmliMuc=meuc(mls  (xindex) )  ; 

%  get  the  data  which  satisfy  x=xrng 
for  isstriiendi 

mx(p)s  mlx(i) ;  mz(p)s  mlz(i);  ms(p)s  mls(i); 

P=p-*-l; 

end 

%  set  the  eucis  of  plot 
IIUUCZ3m2UC(fflZ)  ; 
meuexsmax(xinlroax,xrpomax)  ; 
minx3min(xinlmin,xrpomin)  ; 

2ixis(  [  minx  maxx  0  meucz  ]  ) 
plot (ms,mz, ' -g' ) 
hold  on 

%  plot  the  horizontal  line 
hrs=t] ;hrz=t] ; 
hrs=roinx:  1  imeucx; 
for  izl: length (hrs) 
hrz ( i ) 3anh2 ; 

end 

plot (hrs, hrz) 


46 


%  plot  rpo 

if  xlasl 

plot  (rpsl,  zpzl,  '■►b' ) 

•nd- 

i£  x2332 

plot (zps2, rpz2. * .w' } 

end 

if  x3x33 

plot (  rps3 , rpz3 , 'xb' ) 

end 

if  x43s4 

plot (  rps4 , rpz4 , '*r' ) 

end 

xlabeK'  Propagation  Loss  (dB)  '}; 
ylabeK'  Height  (m)  '); 
gtext(  [nuin2str (prof)  ] ) ; 

gtext  ( [  *< - RPO'  J )  ; 

gtext  ( ( '  < - MLAYER'  ] )  ; 

gtext  ([' range  :  '  .num2str(xrg) , 'Ion' ] )  ; 
gtext  ( ( '  frequency:  ' ,  nuin2str (frq) ,  'MHz'  ] )  ; 
gtext {  ["TX  height:  ' , num2atr (anh) , 'm' ] ) ; 
if  xl==l 

gtext (['FE  region:  +']) 
end 

if  x2==2 

gtext ( ( ' RO  region :  . '  ] ) 
end 

if  x3==3 

gtext (['XO  region:  x']) 
end 

if  x4s34 

gtext  { ( '  PE  region *.  * '  1 ) 
end 
grid 
pause 

csavsinput ( '  save  the  plot  7  (y/n)  ' ,  '  s ' ) ; 
if  (csavss'y' ) I (csav=='Y' ) 
t33('ineta  'rl'.fno]; 
eval(t3) ; 

end 

ans3input('  enter  N  to  Exit  ','s'}; 
if  (ans3='N')|(  euiss='n') 
ansO  ; 

end 

clear  mz  ms; 

clear  rpsl  rps2  rps3  rps4 
clear  rpzl  rpz2  rpz3  rpz4 
clg 

hold  off 

end 
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APPENDIX  B:  PROPAGATION  LOSS  UNDER  THE  INFLUENCE  OF  A  300  M 

SURFACE>BASED  DUCT 

This  Appendix  displays  the  propagation  loss  computed  by  RPO  and  M-Layer  under 
the  influence  of  a  300  m  surface-based  duct  at  3  GHz  and  12  GHz  at  ranges  of  IS,  20, 
30,  40,  50,  60,  70,  80.  90,  100  and  110  km. 

1.  Propagation  loss  at  3  GHz 

Figures  B.l  through  B.l  1  displays  the  propagation  loss  at  3  GHz  computed  by 
RPO  and  M-Layer. 
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Figure  B  J.  Propagation  loss  at  30  km. 


Figure  B.4.  Propagation  loss  at  40  km. 
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Figure  B^.  Propagation  loss  at  SO  km. 
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Figure  B.6.  Propagation  loss  at  60  km. 
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are  B.7.  Propagation  loss  at  70  km. 


Figure  BA  Propagation  loss  at  80  km. 
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Figure  B.ll.  Propagation  loss  at  110  km. 


2.  Propagation  loss  at  12  GHz 

Figures  B.12  through  B.22  displays  the  propagation  loss  at  12  GHz  computed 
by  RPO  and  M-Laycr. 
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Figure  B.14.  Propagation  loss  at  30  km. 
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Figure  B.16.  Propagation  loss  at  SO  km. 


Figure  B.17.  Propagation  loss  at  60  km. 


Figure  B  JO.  Propagation  loss  at  90  km. 


Figure  BJl.  Propagation  loss  at  100  km. 
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Figure  B  J2.  Propagadon  loss  at  1 10  km. 
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APPENDIX  C:  PROPAGATION  LOSS  UNDER  THE  INFLUENCE  OF  A  14  M 

EVAPORATION  DUCT 

This  Appendix  displays  the  propagation  loss  computed  by  RPO  and  M- Layer  under 
the  influence  of  a  14  m  evaporation  duct  at  3  GHz  and  12  GHz  at  ranges  of  15,  20,  30, 
40,  50,  60,  70,  80,  90,  100  and  110  km. 

1.  Propagation  loss  at  3  GHz 

Figures  C.l  through  C.1 1  displays  the  propagation  loss  at  3  GHz  computed  by 
RPO  and  M-Layer. 
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lire  C.l.  Propagation  loss  at  15  km. 
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ire  C  J.  Propagation  loss  at  30  km. 
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Figure  C.4.  Propagation  loss  at  40  km. 
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Figure  CjS,  Propagation  loss  at  60  km. 
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2.  Propagation  loss  at  12  GHz 

Figures  C.12  through  C.22  displays  the  propagation  loss  at  12  GHz  computed 
by  RPO  and  M-Layer. 
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Figure  C.13.  Propagation  loss  at  20  km. 
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re  C.16.  Propagation  loss  at  SO  km. 


Figure  C.17.  Propagation  loss  at  60  km. 
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Figure  C.18.  Propagation  loss  at  70  km. 
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Figure  C.19.  Propagation  loss  at  80  km. 
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APPENDEX  D:  PROPAGATION  LOSS  UNDER  THE  INFLUENCE  OF  A  300  M 


SURFACE-BASED  DUCT  OVER  A  14  M  EVAPORATION  DUCT 
This  Appendix  displays  the  propagation  loss  computed  by  RPO  and  M-Layer  under 
the  influence  of  a  300  m  surface-based  duct  over  a  14  m  evaporation  duct  at  3  GHz  and 
12  GHz  at  ranges  of  15.  20.  30,  40,  50.  60,  70.  80.  90.  100  and  1 10  km. 

1.  Propagation  loss  at  3  GHz 

Figures  D.l  through  D.l  1  displays  the  propagation  loss  at  3  GHz  computed  by 
RPO  and  M-Layer. 
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Figure  Propagation  loss  at  20  km. 
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Figure  D  J.  Propagation  loss  at  30  km. 
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Figure  D.4.  Propagation  loss  at  40  km. 
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Figure  D.6.  Propagation  loss  at  60  km. 
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ire  D.7.  Propagation  loss  at  70  km. 
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Figure  D.8.  Propagation  loss  at  80  km. 
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Figure  D.9.  Propagation  loss  at  90  km. 


Figure  0.10.  Propagation  loss  at  100  km. 
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Figure  D.ll.  Propagation  loss  at  1 10  km. 


2.  Propagation  loss  at  12  GHz 

Figures  D.12  through  D.22  displays  the  propagation  loss  at  12  GHz  computed 
by  RPO  and  M-Laycr. 
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Figure  D.15.  Propagation  loss  at  40  km. 
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Figure  D.16.  Propagation  loss  at  50  km. 
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Figure  D.18.  Propagation  loss  at  70  km. 


Figure  D.19.  Propagation  loss  at  80  km 
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Figure  DJ.2.  Propagation  loss  at  1 10  km. 
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APPENDIX  E:  PROPAGATION  LOSS  UNDER  THE  INFLUENCE  OF  A  150  M 

SURFACE-BASED  DUCT 

This  Appendix  displays  the  propagation  loss  computed  by  RPO  and  M-Layer  under 
the  influence  of  a  150  m  surface-based  duct  at  3  GHz,  6  GHz  and  12  GHz  at  ranges  of 
15,  20,  25,  30,  35,  40,  50,  60,  70,  80,  90,  100  and  110  km.  The  modified  rcfractivity 
profile  is  given  in  Table  E.l. 


Table  E.I.  A  150  m  surface-based  duct 


Z 

meters 

Mi 

0 

0.0 

339.0 

1 

100.0 

350.8 

2 

150.0 

301.3 

3 

1,000 

401.6 

1.  Propagatimi  loss  at  3  GHz 

Figures  E.1  through  E.13  displays  the  propagation  loss  at  3  GHz  computed  by 
RPO  and  M-Layer. 
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Figure  E^.  Propagation  loss  at  35  km. 
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Figure  Ej6.  Propagation  loss  at  40  km. 
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Figure  E.13.  Propagation  loss  at  1 10  km. 


2.  Propagation  loss  at  6  GHz 

Figures  E.  14  through  E.26  displays  the  propagation  loss  at  6  GHz  computed 
by  RPO  and  M-Layer. 
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Figure  E.16.  Propagation  loss  at  25  km. 
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Figure  E.17.  Propagation  loss  at  30  km. 
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re  E.18.  Propagation  loss  at  35  km. 


Figure  E.19.  Propagation  loss  at  40  km. 
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PropagKion  Lon  (dB) 
Figure  EJ3.  Propagation  loss  at  80  Ion. 
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3.  Propagation  loss  at  12  GHz 

Figures  E.27  through  E.39  displays  the  propagation  loss  at  12  GHz  computed 
by  RPO  and  M-Layer. 
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Figure  EJ8.  Propagation  loss  at  20  km. 
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ire  E  J3.  Propagation  loss  at  50  km. 


Figure  EJ4.  Propagation  loss  at  60  km. 


Figure  EJ5.  Propagation  loss  at  70  km. 
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Figure  EJ6.  Propagation  loss  at  80  kiiL 
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Figure  EJ8.  Propagation  loss  at  100  km. 
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Figure  EJ9.  Propagation  loss  at  110  km. 
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APPENDIX  F:  PROPAGATION  LOSS  UNDER  THE  INFLUENCE  OF  A  100  M 

SURFACE-BASED  DUCT 

This  Appendix  displays  the  propagadon  loss  computed  by  RPO  and  M-Layer  under 
the  influence  of  a  100  m  surface-based  duct  at  3  GHz,  6  GHz  and  12  GHz  at  ranges  of 
15,  20,  25,  30,  35,  40,  50,  60,  70,  80,  90,  100  and  110  km.  The  modified  refiractivity 
profile  is  given  in  Table  F.l. 


Table  F.l.  A  100  m  surface-based  duct 


i 

Zi 

meters 

Mi  1 

0 

0.0 

339.0  1 

1 

50.0 

344.9 

1  ^ 

100.0 

295.4 

I  ^ 

1000.0 

401.6 

1.  Propagation  loss  at  3  GHz 

Figures  F.l  through  F.13  displays  the  propagation  loss  at  3  GHz  computed  by 
RPO  and  M-Layer. 
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Figure  F^.  Propagation  loss  at  35  km. 


Figure  F.6.  Propagation  loss  at  40  km. 
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Figure  F.8.  Pn^agation  loss  at  60  km. 
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ire  F.9.  Propagation  loss  at  70  km. 


Figure  F.IO.  Propagation  loss  at  80  km. 
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Figure  F.12.  Propagation  loss  at  100  km. 


2.  Propagation  loss  at  6  GHz 

Figures  F.14  through  F.26  displays  the  propagation  loss  at  6  GHz  computed 
by  RPO  and  M-Laycr. 
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Figure  F.19.  Propagation  loss  at  40  Ion. 
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Figure  F.23.  Propagation  loss  at  80  km. 
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Figure  FJl6.  Propagadon  loss  at  110  km. 


3.  Propagation  loss  at  12  GHz 

Figures  F.27  through  F.39  displays  the  propagation  loss  at  12  GHz  computed 
by  RPO  and  M-Laycr. 
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Figure  FJ2.  Propagation  loss  at  40  km. 
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Figure  F.35.  Propagation  loss  at  70  km. 


Figure  FJ6. 1’ropagation  loss  at  80  km. 


Figure  FJ8.  Propagation  loss  at  100  Ion. 
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Figure  F.39.  Propagation  loss  at  1 10  km. 
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